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Abstract 

We construct a model which describes a recently performed experiment (Phys. Rev. A 64, 
050301(R) (2001)) in which an entangled state between two modes of a single cavity is built. 
Environmental effects are taken into account and the results agree with the experimental findings. 
Moreover the model predicts, for different conditions of the same experiment, a decoherence-free 
subspace. These conditions are analyzed and slightly different experiments suggested in order to 
test its viability. 
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I. INTRODUCTION 



The role of environmental degrees of freedom is a central theme today both in the study of 
Foundations of Quantum Mechanics ll I and for technological improvement in what concerns 

Q 1 

Quantum Computation 2\. Cavity Quantum Electrodynamics |3|] is one of the areas where 
such problem has been intensively investigated in order to establish the dynamics of the 
evolution of quantum superpositions considered as an open system. 

Decoherence control is an important topic in this context. Systems composed by two 
parts interacting with a common reservoir are of interest, since theoretically they may lead 
to the existence of trapped states. This effect shall be associated to cross decay rate terms, 
which have been studied at least since Agarwall j^J. Derivations of master e quat ions for 
such systems may be found in 0, Q, 0, Q and applications of them in jit], [ill I12L \\\ . 

A recent and particularly interesting experiment involves the construction of two electro- 
magnetic field modes in a single cavity In the present contribution we investigate the 
consequences of a straightforward generalization of the Caldeira-Legget type model for the 
discussed experiment and are able to explain the results. Moreover this model predicts the 
existence of a decoherence-free subspace (DFS) provided the coupling to the environment 
satisfies certain conditions. These cannot be tested in the present experiment although with 
slight modifications the feasibility and robustness of such spaces can be easily accessed. 

In Sec. II we present a Hamiltonian to model the system (two electromagnetic field 
modes) subjected to the environment, and derive a master equation in the usual Markov 
regime. A technique for solving this master equation, group theory for superoperators, 
is given in Appendix A. The experiment is described in Sec. Ill, where its results are 
compared to the theoretical ones from our model. In Sec. IV we suggest slight modifications 
in the experimental sketch, useful to investigate a possible tendency of forming DFS in such 
systems. Our conclusions are presented in Sec. V, where the feasibility of DFS is analyzed. 
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II. FROM THE MODEL TO THE MASTER EQUATION 



The Hamiltonian we use to describe two electromagnetic field modes in a single super- 
conducting cavity plus environment is 



where 



H = H + Hint, (1) 

H = H 5 + U E , (2) 
H 5 = M^ajax + M7 2 al,a2, 
He = hy]uJkclc k , 

k 

Hint = ( a lfe a l C fc + a lfc a l C fe) + ("2fca 2 C^ + tt2fc a 2 C A:) ■ 

k k 

The operators ai (a 2 ) and a\ (a|) are annihilation and creation bosonic operators for mode 
M a (Mb) with frequency fli (fi 2 ). The environment is modelled by a set of harmonic oscil- 
lators with creation and annihilation operators cj. and linearly coupled to the system, 
as it is usually done jl^. Harmonic oscillators are appropriate to model extended modes 
. as phonons in the cavity mirrors and electromagnetic environment modes in the lab- 
oratory. The coupling between these environment oscillators and the oscillators of interest 
may occur by complicated processes (e.g., a photon may be scattered from M a to an 
electromagnetic environment mode by an spurious atom inside the cavity). The coupling we 
considered is an effective one, related to changes of one photon between environment and 
system. 

In what follows we deduce a master equation for the system in a similar fashion as done 
in Ref. [jj, where we also considered a coupling between the relevant modes (more detailed 
calculation can be found in Ref. j^l). Let us take p(t) the complete density operator 
concerning the system plus environment. Its time evolution may be given by 

j t p(t) = ^-[n mt (t),p(t)], (3) 

where 

p (t) = ei Uot p (t) e-^ Hot , ei Uot n int (t) e~i Hot 
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are in the interaction picture. The high quality factor of the cavity permit us to consider 
Hint (t) small, since M a and Mf, are weakly coupled to the environment. Disregarding the 
terms of third order in H int (t), Eq. (JHJ) leads to 

~p{t) - p(0) = - l -jjt'[tl int (0 ,P(0)] - Jpf^f dt "^ (0 > [Hint (t"),p(0)]]. (4) 

Let us admit that at t — the system is prepared in the state p$ (0) and the environment 
is in thermal equilibrium. Thus 

p(0) = p s (0)®p E (0), (5) 

with 

Pe(0) = ^II ex P (-^4 c fe) . 2 = [[E ex P^M. /5 = [k B Ty l . (6) 

fe fc n=0 

! iere is Boltzmann constant, T is the absolute temperature and Z is the partition function 
1^ | . Taking the trace over the environment degrees of freedom in both sides of Eq. (jlj), we 
can find, in the limit of zero temperature, 

Ps (t) - Ps (0) = l dt 'l dT | (E «i*«i* e<W * T J e-^ T e i ^- n «)*' (a,p 5 (0) a] - p 5 (0) a 
+h.c, 

where h.c. stands for Hermitian conjugate. 

Notice that J2k a ik a jk eltJkT decays very fast with the growing of r. Thus we may modify 
the integration limits above and obtain 

2 

p s (t) - p s (0) = { + iAij) (aips (0) at - p s (0) aja<) ( f>V^-<^') } + h.c. (8) 

for t ^> r c , where r c is the time within a.ika* k e lu)kT have appreciable values. The constants 
kij and A^- are real, defined by 

hj + iA« = V a ifc a* fc / Ve*^-^. (9) 

Now we differentiate both sides of Eq. (jSJ) and iterate to get, in an analogous way as in 
Ref. H, 

, 2 

J t Ps (*) = E {(*« + ( a *^ (*) a l - & (*) a l a >) + h.c, (10) 



where terms of second order in kij + iAjj, that are of fourth order in Hint (t), were not taken 
into account. Returning to the Schrodinger picture we write the master equation 



where 



j t Ps (t) = C Ps (t) , 



(11) 
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(12) 



is a Liouvillian superoperator (i.e., an operator which acts on operators). We use the 
conventional notation for superoperators [20]: the dot sign (•) indicates the place to be 
occupied by ps (t) , where the superoperator acts. 

In Eq. (I12J) the constants ku and /c 22 are associated to the individual dissipation of the 
modes M a and The constants An and A 2 2 have unitary effects, renormalizing the 
oscillation frequencies fli and Q2 (Lamb shifts). The coefficients A; 12 , &21, A 12 or A 2 i are 
related to a communication channel between the modes mediated by the environment, with 
unitary and non unitary effects over the evolution of the system. In Eq. © we see that these 
terms will be appreciable only if: 1) |aii&| and |a 2 fc| are both not zero for several values of k 
(this means that the system's modes interact effectively with the same reservoir); 2) a\ka\ k 
have phase correlation for different k (the system's modes interact with the environment in 
a microscopic correlated way). Relatively large values of the cross decay terms A42 + i^n 
and k<i\ + iA 2 i are important for the appearance of DFS (see Sec. IV). The experimental 
conditions for it will be discussed in Sec. V. 
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III. COMPARING THEORETICAL PREDICTIONS AND EXPERIMENTAL RE- 
SULTS 



A scheme of the experiment is presented in Fig. [TJ Cavity C supports two modes, M a 
and M b , with orthogonal polarizations and different frequencies. The M a mode frequency 
is larger than that of the M b by 5. Circular Rydberg atoms A s and A p , with levels called 
e and g, are sent through the cavity, the first one to create the entangled state and the 
second one to reveal its quantum nature. The detuning A between the e — ► g transition 
frequency and the mode M a may be adjusted by means of the Stark effect. The sequence 
of the experiment is described below, closely following Ref. |l4j, where all environmental 
effects have been ignored. 

At time t — 0, the atom A s enters the cavity in state |e). The state of the atom- two 
modes system is then |e, a ,0&). The parameter A is initially set to zero and the atom 
interacts with M a through a | Rabi pulse. The interaction with M b is ignored due to the 
detuning and the atom-field coupling is considered to be a constant (vacuum Rabi frequency 
Q) . With an appropriate phase choice of the atomic dipole and assuming the energy of the 
state \g, l a , b ) as zero, the atom-cavity state at t = ^ will be 



* It 



-) 



V2 



[|e s ,0 a > + \9sAa)] |0 6 > 



The atom is next, by Stark effect, in resonance with mode M b (A = —5). If one now neglects 

id tal 

3ttY 



the interaction with mode M a and takes care of the phases appropriately, one gets for t = |^ 



(13) 



where (ft = f + ff • The state (|g)) of the atom A s ends up factorized and need not be 
considered any longer. 

At time T (T > |^), atom A p enters the cavity, whose state will be given by 



i» ( *=r» = ^ 



i exp 



exp (iST) \0 a ,l b ) + \l a ,0 b ) 



A p is initially in the ground state and interacts through a 7r Rabi pulse with M a . Next, it 
interacts with mode M b through a | Rabi pulse, yielding the state 



i5T e i5-K/2Q. 



)]■ 



Mode M a ends up in the vacuum state, factorized. The probability P e {T) of finding A p in 
state |e) is therefore given by 

[1 + cos (8T + $)] 



(14) 



where 



Sit 
2XV 



In Appendix A we present a general solution of the Master Equation (fTT| . This may 
be used to describe the experiment in question taking dissipation into consideration in the 
period between the crossing of the two atoms. In fact, assuming that in t — |S the cavity 
state is given by Eq. ()13|). in t = T, considering the interaction with the environment, will 
be 



Ps(T) = -[(e^F 2 (r) + L 2 (r))|0)|l) + (F 1 (r) + e^L 1 (r))|l)|0)] [h. c] 



+ 1- 



e^F 2 (r) + L 2 (t) 



\Fi (r) + e^Li (r) 



|0> |0) <0| (0| 



where r = T — f£ and 
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(15) 
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r = ^c 2 + (fci2 - iAi 2 ) (A; 21 -iA 21 ), 



1? 



fcn + /c 22 , (O x - A n) + (0 2 - A 22 ) 

— r~ ^ 



2 2 
Using the above solution for p s {T) it is a simple matter to evaluate P e (T). We get 

P e (T) = i |- (iq (r) + e^Lx (r)) + *e 2t$ (e*F 2 (r) + L 2 (r)) f . 

The maximum possible values for \k± 2 \, \k 2 i\, |Ai 2 | and |A 2 i| are of the order of ku and k 22 . 
Given the large detuning, the experiment is not sensitive to these cross decay constants and 



is consistent with k±2 = &21 = A12 = A21 = 0. The experimental results are also consistent 
with An = A 2 2 = 0. In this case we get 



2k\\T _|_ e -2fc22T 



+ e -(fcll+fe 22 )r cos (^ T + $) 



(16) 
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Making k\\ = k 2 2 = 0, we recover expression (frlj) . 

Note that Eq. has been derived for zero temperature. This is not the exact experi- 
mental condition since in thermal equilibrium the modes M a and M b contain a small fraction 
of thermal photons (~1). In order to take this effect into account we use effective dissipation 
constants for the modes k\{ = and k%, = where n is the average number of thermal 
photons and T r a (T r)b ) the measured decay time for mode M a (M b ). 

There are several sources of imperfection in the experiment, related to the construction of 
the modes and to the atomic detection. The fidelity of the state in Eq. fllHj) is estimated to 
be of 50%. Therefore the interaction with the environment is not the only source of visibility 
loss. In order to compare our model results to the experiment we use a reduction factor of 
50%. 

In Fig. El we show P e {T) x T for the same time windows as in Ref. j 141 ] . Eq. (fTB*|) 

have been used and the experimental reduction factor was taken into account. There is 
good agreement between our predictions and experiment, especially in what concerns the 
amplitude and period of oscillation. We note a slight shift, approximately constant in all 
cases (< between theory and experiment. Some possible sources of this phase shift are: 
1) The assumption of a constant atom-field interaction. 2) During the time of switching the 
atoms off resonance with mode M a and in resonance with mode M b , the phase accumulation 
happens in a way that depends on the details of the process, not included in this model. 
Notice that these times are, accordingly to Q|, within 1/zs, ie, in the same order of the 
shift. 3) We have assumed that the atom interacts with just one mode at a time. The 
simultaneous interaction of the atom with both modes create a communication channel 
between the modes, which has consequences on the phase in P e (T). This effect may be more 
relevant during the switching of the atom, when the frequency of the atomic transition is 
not maximally far from the frequencies of M a and M b . 
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IV. SEARCHING FOR DECOHERENCE-FREE SUBSPACES 



We next perform a mathematical analysis of a different situation than the experimental 
one: let us consider resonating modes. In this case 5 = and the role of A; 12 and fc 2 i 
become significant. In Fig. El we plot P e (T) x T for the case, 5 = 0, k 12 = &21 for various 
values of k\2, keeping the other parameters as considered in Sec. III. Notice in Fig. El that if 
ki2 — \Jk\\ki2 the field in C does not go to zero for long times, suggesting a decoherence-free 
situation. In fact in this situation the constant A_ (Eq. (|15j)) has a vanishing real part and 
the exponentials related to it will therefore not decay. The field decay is then completely 
dictated by A + . So part of the field remains protected. 

What does this condition /c 12 = \fk~\\h^ mean physically? It means that the environment 
acts as a coherence feed-back mechanism. It means that photons would scatter from one 
mode and be transferred to the other without loosing their coherence. This may seem rather 
unrealistic although it is a sound mathematical consequence of the extension of a model 
which works very well for a single mode Moreover, the other curves in Fig. El indicate 
that even very slight deviations from this condition already destroy the existence of this 
decoherence-free situation. 

The experimental scheme to investigate the possibility of such an effect would need one 
as small as possible detuning. Of course the assembly must be altered if we don't want the 
atom interacting simultaneously with both modes. We have at least two strategies: 

1) In [14| a quadratic Stark effect is obtained by applying a dc voltage across the mirrors, 
which maintains the atomic orbital plane perpendicular to the cavity axis. In this case the 
atoms couple equally to both modes. We may get the Stark effect by, instead, applying 
two dc voltages perpendicular to the cavity axis (and perpendicular to each other). Thus 
the atomic orbital plane may be maintained perpendicular to the M a or polarizations, 
and the atom may be coupled to just one mode at a time. Such voltages may be produced 
directly in the ring around the cavity (in this case the ring must not be continuous) or, if 
we take away the ring, in plates outside the cavity. 

2) The modes may be constructed in two separate cavities. Since we are investigating 
effects of the interaction between the cavities modes through the environment, it would 
be suitable to have the cavities as close as possible to each other. In fact, if the distance 
between the cavities is small in the modes' wavelength scale, this interaction is expected to 
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be maximized. 

Adjusting the curve P e {T) could teach us something about the values of k\ 2 and k 2 i, or 
at least about a tendency to form a DFS, if one approaches the ideal limit described above. 

From the theoretical point of view, it is an important issue to be able to select a specific 
DFS. Consider that M a and Mf, have the same frequency u) and 

k 22 + iA 22 = K 2 (k n + iA u ), (17) 
k 12 + iA 12 = k 21 + iA 21 = k (ku + iA n ) , 

where k is real (a particular choice of k is related to the quotient between the quality factors 
of the modes). For all Ay = this is the case treated in this section when k\ 2 = k 2 i = 
\Jk\\k 22 . Defining the bosonic operators 

A = 1 (ai + «a 2 ) , 

V 1 + K 

B = (a 2 - /cai) , 

V 1 4" K 

we write the Liouvillian (|12|1 as 
C = Ca + ^b, 

^ = ~\ (Pa, •]) + (! + (2A . At - A+A . - • A 1 " A) , 

>Cb = 4([H B! .]), 

H A = h(u — A aa (l + k 2 )) A 1 " A, 
H B = ^B f B. 

Thus, if a the system is in a state which may be written as 

Ps = 5> m ,„ (BT |0, 0) (0, 0| B m + h.c, (18) 

m,n 

it is not affected by the environment. The states (|18|) define a DFS. Relevant examples for 
Quantum Optics and Quantum Information are the coherent state 

Ps = \-kv) \v) (h.c.) , 

the superposition of coherent states 

p s = N (\-kv) \v) + e**\-K,w) \w)) (h. c.) 
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and the superposition of Fock states 

PS = ^=L=(|O a >|l b )-K|l a > |0 b ». 

v l + « 

We shall emphasize that since the conditions (|T7|) are related to the characteristics of the 
environment, they can not be freely chosen by the experimenter. To understand physically 
what may lead to conditions (|17|). consider that the coupling constants of M a and Mb to the 
reservatory modes may be factorized in the form 

= aiTfc, (19) 

where «i and «2 are real numbers. This corresponds to M a and Mj, interacting with the 
environment in a microscopic correlated way, with a possible difference in the intensity of 
the interaction |2lJ. For ressonant modes, conditions (fTTIJ) imply conditions (jl7j) (see Eq. 
(JIU)), and we get the DFS just described. Another way to relate conditions to this DFS 
is to use them directly in the Hamiltonian (j2J) 2^|. An experimental sketch to lead to this 
microscopic correlation would need modes as close as possible to each other, preferentialy 
with the same polarization, as will be discussed in more detail in the next section. If the 
correlation achievable in an experiment is not perfect, ki2 + zAi2 and foi + 1A21 will assume 
intermediate values between zero and the ones in (fT7|) . 



V. CONCLUSIONS 



We deduced a master equation for two oscillators in the presence of a common reservoir 
and used it to model a Cavity Quantum Electrodynamics experiment involving two modes 
constructed in the same cavity. The theoretical results show good agreement with the ex- 
periment. Such a master equation predicts the existence of DFS if its cross terms (terms 
involving operators of both oscillators) have sufficiently large coefficients. Since the exper- 
iment analyzed is not sensitive to these coefficients, we proposed two slightly modificated 
experiments which may permit to investigate them. 

The model indicates that a DFS may appear if M a and M& modes interact with the 
environmental modes in a microscopically correlated way. Of course it will not be the 
case for most systems, and it is not a simple situation to construct. Probably we may 
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have this microscopic correlation at least partially for modes whose distance in space is 
small compared to the wavelength of the most important environmental modes (the ones 
with frequencies near the frequencies of M a and M b , as may be seen in Eq. Q). If the 
environment is composed mainly by electromagnetic modes, the relevant scale is the scale 
of the wavelengths of M n and M&. 

The experiment in was performed with modes in the same place in space. Unfortu- 
nately it doesn't guarantee the microscopic correlation, since M a and M& have orthogonal 
polarization, and then they "perceive" different microscopic environments. Although it must 
be not easy to built, there is no theoretical impossibility to construct modes close in the 
scale of their proper wavelengths. 

The model states clearly that it is very difficult to achieve the parameters necessary to 
observe a DFS, and Fig. |3] shows how fast the DFS is spoilt when we leave the perfect 
situation. It is a sign that this model is realistic. But the model also indicates where are the 
main difficulties and what may be done to approach the ideal conditions. Due to the Cavity 
Quantum Electrodynamics current experimental stage, the sketches we propose may just 
investigate the tendency of forming DFS. If this tendency is confirmed, it shall encourage 
later developments. 
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APPENDIX A 

Although thoroughly derived in Ref. j^, we repeat here the solution to the master 
equation, with slight modifications, for completeness. We will use the parameter derivation 
technique, which allows one to determine coefficients q such that the identity 

e (7iOi+720 2 +-+ 7n O„)t _ e a(t)Oi e ?2(t)0 2 . . . e <r„(i)O n ^ 1) 

is valid, where the Oj's are superoperators forming a closed Lie algebra and t is a parameter. 
The parameter derivation technique consists of the following procedure: 
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1. Derive both sides of Eq. (|A1|) with respect to t and get 

' n \ / n \ n 

^ 7i O, exp ^ 7i CM = & (t) Ox H e*W°< 



i=i / \j=i / i=i 



+ft(t) e «(*)°i0 2 jie«W°' 

i=2 

+ ••• 

n-1 

(t) IJ e ? ' (<)0 'O n e ? " wo ". (A2) 

i=l 



2. Use the similarity transformation 

<'° G;< '° = e B ^'* ] 0< (A3) 

and the linear independence of {Oj} in order to obtain differential equations for the 
parameters ft (£). 

3. Solve the c- numbers differential equations and obtain the factorized evolution super- 
operator written in the right hand side of Eq. (jAl|) . 



As an example of step EJ above, let us take the second term in the r.h.s. of Eq. (jA2|) : 

e a(*)°i0 2 = e ?l{t)Ol 02e" ?l{t)Ol e tl( * )O1 
= e ?(t)[Oi,«]Q a(*)Oi 

ee h (ft (f) , {O,}) e*(*><\ (A4) 
Analogously one can carry out a similar operation for all the other terms, define 

/ 3 = / 3 (ft (*) , ft (t) , {O,}) , • • • , f n = f n (ft (t) , ft (i) , • • • , ft_x (t) , {O,}) (A5) 
and write 

^7iOi ) exp (^TiOit j =(?iOi+? 2 / 2 + .--^/ n )exp(^7 < 4 t) . (A6) 

V i=l / \i=l / \j=l / 

Equivalently: 

n 

7*0 4 = ftd + ft/ 2 + • • • Un- (A7) 



i=l 
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Due to the linear independence of {Oj}, one obtains a system of coupled differential equa- 
tions for the q's comparing the coefficients of each Oj. 
In our case we have 



Ps (t) = e ct Ps (0) 

_ e fti(t)ai«a{ e /i2(t)a2»a| e z;(t)ai»a| e z(t)a 2 «a{ e ni(t)«a{a2 e n(t)a|ai« 

e m2(t)a^a2» e P2(t)»a^a2 e mi(t)aJai» e pi(t)»a{ai e g(t)a{a2» e g;(*)» a 2 a ip 5 , (Q) (A8) 

Using the method just described we get 

mi (t)-n(t)q(t) e mi «- m2 M, 
,- 22 = m 2 (i)+n(t)?(t)e mi W- ma (" 1 

n (t) + n (t) (mi (*) - m 2 (*)) - n {tf q (t) e mi{t) ~ m2{t \ 
hi = i>i(t)-n l (t)q l (t)e? 1 ®-'»®, 
p 2 (t)+n, (i)^*)^"^, 

c 2 i = (t) + m {t) (p! (t) - p 2 (t)) - n ; (t) 2 g, (t) e Pl «- p2 « 
en = z (t) (iA 21 - k 21 ) 

-zi (t) (iA 21 + k 21 ) + /i! (t) (-2fc n ) + h\ (t) , 
2k 22 = zi (t) (iA 12 - k 12 ) 

-z (t) (iA 12 + k 12 ) + h 2 (t) (-2k 22 ) + h 2 (t) , 
i (A 21 - A 12 ) + k 21 + k 12 = z (t) (i (Q x - A n - fi 2 + A 22 ) - k n - k 22 ) 

-h 2 it) (iA 21 + k 21 ) + h (t) (iAia - k 12 ) +z(t), 
i (A12 - A 21 ) + fc 12 + fc 2 i = (t) (i (fi 2 - A 22 - fii + A n ) - A; 22 - fc n ) 

-/ix (*) (iA 12 + k 12 ) + h 2 (t) (iA 21 - k 2l ) + zi it) . (A9) 



i (A n 






i (A 22 - 


-n 21 ) 


-A; : 




tA 12 


-it 




iA 2 i 


- h 




-A n ) 


-k 


i{n 2 - 


- A 22 ) 


-h 




-lAia 


-k 




-iA 21 


-h 






2k 



12 
21 



22 
12 
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The solution reads 



e mi(t) = Fi ^ ; e m 2 (*) = e -2Ht e -mi(t) } 

M*) = (|F 2 (t)| 2 + |L 2 (t)| 2 )e 4fc -'-l, 

h 2 (t) = {iF^tf + lL^t^e^-l, 

z(t) = - {L x (t) F* 2 {t) + L* [t) Fx (t)) e 4fc -', 

*i (*) = z* (t) , 

n,(t) = (n (t))*, ®(*) = (g (*))*, 

p 2 (t) = (ma (t))* , pi (t) = (mi (*))* , (A10) 



where Fi (t), F 2 (t), L\ (t) and L 2 (t) are given by Eqs. (fTHj) . 
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FIG. 2: Probability P e (T) for detecting A p in state e for the experimental situation described in 
Ref. [ljj. The curves refer to the dissipative model. The dots are from the experiment. 
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FIG. 3: Probability P e (T) for detecting A p in state e for resonant modes, according to the dissi- 
pative model. We consider several values for k±2 = k2i- ki2 = hi\ = 0, k\2 = k2\ = 0.5^/^17^22, 
ki2 = ^21 = 0.7 v / fcii^22, ku = k2i = 0.9y / ¥nlc22, k\2 = &2i = ^Aii^22- Upper curves correspond 
to higher values for ku = &2i. 
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